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Abstract 



The overdamped dynamics of a charged particle driven by an uniform electric 
field through a random sequence of scatterers in one dimension is investigated. 
Analytic expressions of the mean velocity and of the velocity power spectrum 
are presented. These show that above a threshold value of the field normal 
diffusion is superimposed to ballistic motion. The diffusion constant can be 
given explicitly. At the threshold field the transition between conduction and 
localization is accompanied by an anomalous diffusion. Our results exemplify 
that, even in the absence of time-dependent stochastic forces, a purely me- 
chanical model equipped with a quenched disorder can exhibit normal as well 
as anomalous diffusion, the latter emerging as a critical property. 

PACS number(s): 05.40.-a 



I. INTRODUCTION 



Various dynamical models have been introduced for describing anomalous diffusion. A 
good deal of them rely upon a stochastic formulation of the dynamical rule, as in the 
celebrated paper by Ya. G. Sinai, describing a random walk in a Id random potential [|TJ. 
More generally, the dynamics is either given by a master equation, with random hopping 
rates, describing the effect of a random environment, or, in the continuous version, by a 
Langevin equation with a time-dependent random force, usually assumed as Gaussian. The 
literature in this field is so huge that here we limit ourselves to mention the excellent review 
article by Bouchaud and Georges 0. 

Anomalous as well as normal diffusion have been tackled also in the framework of deter- 
ministic dynamics. A well-known example is given by the standard map ]3H, where normal 
diffusion is observed when the phase space is represented by a chaotic sea with sparse sta- 
bility island QJ^]. Conversely, when many stability islands coexist finite time trapping of 
chaotic orbits around them induces correlations yielding anomalous diffusion ||. Strong 
anomalous diffusion has been observed in other deterministic systems: here we just mention 
two models of Id intermittent maps introduced by Geisel et al. and by Pikovsky |7|]. In 
the former example the complex scenario of normal and anomalous diffusion emerging from 
a chaotic dynamics was pointed out for the first time. In the latter example, it has been 
shown that one can pass from normal to anomalous diffusion while varying the polynomial 
behaviour at the unstable fixed point of the intermittent map. It is worth observing that all 
the above-mentioned examples of deterministic dynamics exhibiting anomalous as well as 
normal diffusion concern chaotic maps, i.e. time-discrete unpredictable dynamics. Another 
example of diffusion created by a deterministic dynamics, that of a one-dimensional Brown- 
ian particle subject to elastic collisions with 'light' particles was given in ||. Depending on 
the asymptotic scaling of the mass ratio between the test particle and the other particles, 
the asymptotic process was shown to be either Ornstein-Uhlenbeck or Wiener. 

In this paper we consider the purely mechanical problem of an overdamped charged 
particle moving along a line and submitted to an electric field and to a random potential 
created by a set of quenched random scatterers. At variance with similar models that have 
been recently investigated P,[10| , we do not include any random time-dependent force of 
the Langevin type. 

We show that if the electric field is larger than a critical value, a current is created and the 
particle exhibits, on top of the ballistic motion, a diffusive behavior. The diffusion constant 
can be explicitly computed. It is a function of the average distance between scatterers and 
of the mean value and the variance of the passage time through a scatterer, and it vanishes, 
as expected, with vanishing randomness. Below the critical value of the electric field also the 
current vanishes. The transition from an insulator to a conductor induced by the electric 
field can be described by analogy with phase transitions as being first or second order. In 
the former case, the current and the diffusion constant do not vanish and remain finite at 
the threshold, whereas in the latter case the current vanishes and the diffusion constant 
diverges according to a power-law behavior. The critical properties can be characterized 
by a scaling law for the power spectrum of the current. Our study shows that quenched 
disorder is sufficient to create normal diffusion and also anomalous diffusion, at least for the 
critical value of the electric field. 
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We note that our model can be interpreted as the one-dimensional version of the following 
problem: a massive particle submitted to a strong friction is sitting or sliding on a rough 
surface, inclined with a certain angle a with respect to the ground. This model problem 
seems to be important in understanding the phenomenon of segregation by flow. A variant 
of it was recently studied numerically and through a stochastic model | IT| . In our treatment 



the component of the gravity force parallel to the surface plays the role of the driving field. 
The irregularities of the surface, as felt by the particle, are modelled by a random potential. 
It is only when a passes through a certain critical value that the particle will be able to fall 
down indefinitely, and then its motion will exhibit a diffusive behaviour, that turns out to 
be anomalous at the critical value of the angle. 

In the next section we derive the general formulae which will be used in section 3 for 
solving two particular models. In the first model, that we call a gas (or, better, a glass, 
since frozen in), the scatterers are distributed uniformly in an interval. In the second model, 
which can be considered as a special case of the first one, and called a crystal, they follow 
each other at equal distances. In section 4 we discuss the transition between conduction and 
localization through a specific example. The derivation of a scaling function for the velocity 
power spectrum and some concluding remarks are contained in section 5. 

II. GENERAL SETUP 

The equation of motion of an overdamped particle in a constant external field E > is 

x = E - V'(x) (1) 
Throughout the paper we suppose the potential to be a superposition of disjoint scatterers, 

N 

V(x) = (f>j(x — Tj) , 4>j(x) = if |x| > a (2) 
j=i 

with centers rj in an interval [L ,Li], subject to the constraints 

rj + i — Tj > a > 2a (3) 
This also implies L\ — L > (N — \)a. We introduce the definitions 

rf = Tj ± a (4) 



and 



and suppose 



tf : x(tf)=rf (5) 



tj < tf < tj +1 (6) 

that is, the particle moves steadily from left to right, 'above' the barriers. This imposes a 
condition on the external field, namely 
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E > 4>'j(x) for all j and x (7) 

The particle starts at t = in x(0) = r < rf and for t < t7 it moves with a constant speed 
x = E, so that 

= r + £t (8) 

By solving for t we obtain, 

If t7 < t < fj , i.e. r~ < x < Tj~, the equation of motion reads 



dt 

whose solution is 



-(x-r j )=E-<f>' j (x-r j ) (10) 



x—rj 



In particular, the passage time above the j-th scatterer is 

=*'<»> = £*4oij (12) 

The solution between successive bumps rt < x < r7 +1 , i.e. <t < tj +1 again corresponds 
to the free case x = E, so that 

t- +l - 1* = Sti^i da) 

Finally, for t > we are again in the free case and integration yields 

x = r+ + E{t - t+ ) (14) 

For the running solutions that we are looking for, the overall equation of motion can be put 
in the form 

x = E-Y J - ri)x (t - *+>(*) (15) 

where x denotes the characteristic function. The Laplace transform of this equation reads 

roo TP N „t+ 

/ x(t)e-i*dt = - - V / J e-ffiAxti) - rMt (16) 
Jo \i j=1 Jtr 

Substituting <p'Ax — Tj) from ([TOD an d using the definitions (|TT|), (0) an d 



^•(/i) = r e~>*'Wdy (17) 



we obtain 



or, with 



E 



A' 



X 



or,- -(! - /; " 



C^/i) = /ic^/i) - E(l - e~^) 



(18) 



(19) 



H / x(t)e"^dt = E + V C, (/x) 



(20) 



Time average of the velocity will be obtained by sending first iV to infinity and then fi to 
zero in the last equation. Since the distance freely run over by the particle up to time tj is 
Tj — r — a — (J — l)2a, we obtain 



Tj (2j - I) a + r 



k=l 



E 



E 



(21) 



We can take Tj random or not: if they are random the condition (|3|) must be fulfilled. Also 
4>j can be random or not: if they are random we suppose they are identically distributed and 
independent, and also independent of all rv By averaging ( |20"D over disorder and noticing 
that e - ^ is independent of Cj(fi), and Cj(fx) are identically distributed, we obtain 



N 



H \ x(t)e-^dt = E + CU) Y 
Jo j=1 



E + B N {n) 



(22) 



In what follows, we drop the subscript of averaged quantities whenever the average is inde- 
pendent of the subscript. By using (^) we can write 



(23) 



The average factorizes and becomes a power because and hence e~ flTk are independent 
and identically distributed. 

A relevant quantity we are looking for is the time-autocorrelation function of the velocity. 
The average of the Laplace transform of this quantity can be inferred from equations (|20|) 



and (|22|) and reads 



fifi e 
>o Jo 

N 



-fit-fi'f 



x(t)x(t') -x(t) x(t')]dtdt' 



£ e-^e~^'Cj(fi)C f (^') - C{ji)e~^' C(jj,') 

N 

CV)CV) J2 + A(//, //) + A(//, fi) - B n (ji)B n (ji') 

i=i 



(24) 
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Here 



l<j<j'<N 

= C(j7) e-^ T C(/i) e^ /E)l{2j ~ 1)a+ro] e^' /E)l{2j ' ~ 1)a+ro] 

1<3<3'<N 



e -(»Wr r )/E D ^ + ^y-i D^y'-j-i (25) 



III. MODELS 

A. A gas of scatterers 

In the first model we are going to consider the scatterers are uniformly distributed in 
the interval [L ,Li], but respect the constraints @. We introduce a new set of random 
variables {xj} through 

rj =L + (j - l)a + Xj (26) 
The joint probability density of x±, . . . , xn is chosen to be 

N\ N 

p{xi, . . . , x N ) = — ^ Y[ 6{x j+ i - xj) (27) 

L 3=0 

where xq = and xn+i = L := L\ — L — (N — l)a. Then p is two- valued and nonzero if 
and only if < x\ < X2 < • • • < xn < L. 

To compute the average over xi, . . . , xn it is useful to introduce 

Pj (x) := 6(x - Xj ) = 7^ (— T)j V^J)T (28) 



and, for j < j' 



N\ x j ~ l (x' - xy'-'- 1 (L - x') N ~ j ' 



Pjf (x,x') := S(x - Xi)S{xf - x r ) = — j— ^ {N _ f)] (29) 



N ,L 

J2 e-**i = e~» c/E J o dx e~» x/E ]T Pj {x)e- p ^-i) ( 30 ) 



First we calculate the mean asymptotic velocity. From (|23|) and (^6]) we obtain 

N rL 

3=1 

where c = L — a — r and 

e -/3(M) = e -^(^-2«)£)(^) (31) 

We can choose r = L — a, and thus c = 0, without restricting generality Summation over 
j can be performed by the use of the binomial formula. It yields 
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A' 



e ^ ' 



f dxe-^[l-^(l-e- /3 )l 
Jo L 



(32) 



We send iV and Li — Lo to infinity so that the mean distance £ exists and I = lim(Li — 
L )/N > a. This yields lim N/L = (£ - a)" 1 and 



E 

With (E2I) and the notation 



(£ - a) + 1 - e 



(33) 



B(n) = \imB N (n) = C(ji) 



{£-a) + l- 



we can write 



roo 

H \ x(t)e^ 1 = E + B(fi) 
Jo 

When fi goes to zero we asymptotically find 



(34) 



(35) 



C(fi) = fi[2a - Er] 
P(H) = H[T+ (a-2a)/E] 



and the limit of (|35D is therefore 



Et 



x 00 



£ + Et- 2a 

This is the time average of the velocity over an infinite run, that is, 



1 r T 



x(oo) = \im — I x(t)dt = Jim 



x{T) - x(0) 



Remarkably, this is different from the average of x(x) over an infinite distance, 

1 rrg+r 



v(x) = lim - T° \e - V'(x))dx = E 

r^oo r J 



(36) 



(37) 



(38) 



(39) 



which is the same as in the absence of the random potential. The reason is that the work 
of each scatterer exerted on the particle is vanishing, 

<j> j (r+ + 0)-<j> j (rJ-0) = 0. 



Notice that the result fl37|) could have been obtained without any computation: x(oo) is the 
average velocity over any interval of length I containing (the support of) a single scatterer. 

To compute the time-autocorrelation function of the velocity we need to evaluate 
A(/i,/i'). We first rewrite it as 



A(/i, fi') = C(n') e-^C(fi) Yl e-&°-*°W-%-'k(*-2 a )<J'-i) 
( L 'dx [ L dx' Pjj ,(x,x')e-E x e ~^ x ' Dbi + fi'y- 1 

JO J x 



(40) 



With fl3l 



A(/i,/i') = C{p') e-»' T C(fi)e~^ (a ~ 2a) dx dx' 



Pjf(x,x')e 

i<j<j'<N 



e E e E 



-/31jrt.(S)(j-l) e -pQ/)V'-j-l) 



Inserting pjjt from (|29|) and using the trinomial formula 

Ml 



E 



li.m—1 „M—m 



a l b m - l c 



(a + b + c 



M 



<i< m <MlKm-l)\(M-m)\ 
with Z = j - 1, m = f - 2, M = iV - 2, 

a = e -/3(/^') x ) 5 = e -/V) (x'-x), c = L-x' 

we get 



A(/i,/i') = 

L r L 



L 2 



C{p') e~^C{p)e 



X 



dx / dx'e ^ e ^ (l + — 
o Jx I L 



In the limit of N and L going to infinity this yields 



x + e'^Xx' -x)-x' 



N-2 



A(p,p') = {Z-°Y 2 C(m') e-^C(/x)e-^ (CT - 2a) 

/■°° _m _jj , ( 1 

dx / dxe eX eX exp 

o Jx 



a 



e -f3(ji+/j,') x _|_ e -/3(/0|V — 



x — x — X 



Performing the integral we finally arrive at 

C{p7) e~^C(p) e-E(-^) 



A(p,p') 



%{e-<r) + l- - a) + 1 - e~^+^ 



e-»' T C(p)/C{p + n') \ B(fjf)B(fj, + p!) 
For two random processes f(t) and g(t) we introduce the notation 

/•oo /-oo 

K f]9 fayt) = ] j Q dtdt'e-^-^ifitM^-f^git')] 



(41) 



(42) 



(43) 



(44) 



(45) 



(46) 



provided the double integral exists. From equation (^4j) and the subsequent computation 
we find in the limit of N going to infinity 



Ki\i(fi, p') 



B(p + p') 



pp' C{p + p'] 



C(p)C(p') + e -^ {a - 2a) e-»' T C{p)B(p') 
B(p)B(p') 



Mi 



(47) 



Set now p = e/2 + iu and p' — pf. The velocity power spectrum is 



+ 2Re 




e-^[x{t) -x(t)]dt 



C{iu) ]} 



Cj-t^e^-^Cjtu) 
_ e i%(°-A) _ _ a) 



where 



A = 2a - Et 



(48) 



(49) 



We recall that C(/i) was defined through equations (pjj), fll7|) and ([TT]) , and in the last 
line of (HS[) the bar stands for averaging over the remaining (single-scatterer) randomness. 
Equation (f48|) is the main result of our paper. For u real, S±\±(u) is a real, even and 
nonnegative function. If E > E c = sup Jj3 , (f>'j{x), the passage times Tj are distributed on a 
bounded support, the velocity correlations decay exponentially and S ± \ ± (u>) is also analytic at 
to = 0. In this case the asymptotic displacement of the particle is a drift with a superimposed 
normal diffusion. Indeed, the diffusion constant D is given by 



2D 



lime^ 

40 



x{t) - x(t)] 2 dt 



lim 

40 Jo 



eK 



x I X 



(f(l+^),|(l-zy))dy 



i + y 2 



71 



S±\±(Q) 



(50) 



Notice that eK^ in the integrand has a uniform upper bound. By using equations 
(p3|), (pl|), (|33| ) and (|36D, a somewhat tedious computation yields 



2D = 5^(0) 



E 



(£ - A) 3 



e 2 (A 2 -A 2 ) + (£-a) 2 A 2 



(51) 



Thus, D depends on the randomness through the averages £, r and r 2 . We conclude that 
in our example of a deterministic dynamics, with a uniformly distributed set of quenched 



random scatterers, there is normal diffusion. If E 
diffusion. We will discuss this phenomenon later on. 



sup 



. x 



we can loose normal 



B. A crystal of scatterers 



Here we choose rj = L + (j — l)er, that is, the scatterers are placed equidistantly. Because 



now rj + i — rj = a, in the limit when N goes to infinity we also obtain 
directly be substituted in 



, (H) and (|5D to obtain 

xioo) = 

v ; a + Er-2a 



a. This can 



(52) 



E 



a — A u 



= — { |CM| 2 + 2Re 



C(-iu) eW'-^Cjiu) 



(53) 
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and 



A 2 — A 2 

5,1,(0) = Ea^j—^ (54) 

We note that substitution of t = a in equations ( pT7| ) and ( |4"8"[ ) eliminates, at least asymp- 
totically, the randomness of (still admitting xn = o(L)). If we also drop averaging over 
the potentials <f>j, all randomness is lifted, and by the general definition fl4"S|), K ± \ ± and S±\± 
must identically vanish. It can be easily verified that formulae (|47|) and (f48|) or fl53|), and 



also (51), indeed show this property. 



IV. FROM CONDUCTION TO LOCALIZATION 

When E varies continuously and passes the value E c = sup JiX (j)'j(x), it switches between 
a conducting state for E > E c and an isolating, or localizing, state for E < E c . At E = E c 
there may be localization, if r = oo (see equation fl37|) ) , and there may also be conduction 
if r remains finite. In both cases there may still occur a large variety of different situations, 
characterized by different critical exponents and normal or anomalous diffusion. The tran- 
sition between conduction and localization bears a resemblence with phase transitions. For 
instance, E can be considered as the analog of the temperature T, the region E > E c that of 
T < T c , and x(oo) may correspond, e.g., to the spontaneous magnetization and the diffusion 
constant to the static zero-field magnetic susceptibility. We may have first and second-order 
transitions and varying exponents depending on the form of probability distribution of the 
passage time r. 

Let us discuss a simple example which illustrates the different possibilities. We consider 
4>j(x) = fj x (|a;| — a) , < a (55) 

so that 4>'j(x) = ifj- Then 



1 



Tj = a 



+ -ir-r (56) 



E-fj E + f. 



Suppose that the common probability density of the random forces fj has a bounded support 
[b, c], where c > and |6| < c. Boundedness is needed to have a transition, and with the 
above choice E c = c. Let us consider on this support a one-parameter family of probability 
densities 



*(«> - < ^^ 1 <-) 



(c-6)- 



with 7 > — 1. For E > c, 



b^ +1 Jo 



c-b 

V 1 



1 



E + c — v v + e 



n 

dv (58) 



where we have introduced the short-hand notation e = E — c. As e goes to zero, asymptot- 
ically 
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(7 + IK r- 





■ i r 


/ V 1 




/o 


2c — v v. 



dv + o(l) , n < 7 + 1 



7+1 



(7 + i)a 

(c - ; 
(7 + l)a n 



c - 6 

In h 0(1) , n = 7 + 1 



(c-6) 



7+1 



n — 7 



+ 0(e 



-n+7+2\ 



n > 7 + 1 



(59) 



In particular, 



(7 + I) a f c ~ b 



(c - Jo 



1 1 

2c — v v 



■^-rln— + 0(1) , 7 = 



dv + o(l) , 7 > 



|7|)« 



and 



|7|(c-6)H7 
(7 + l)a 2 



171 +0(e 1 - hl ) , 7 < 



(60) 



(c - b)^ 1 Jo 

(2a 2 
(c-6)= 
(7 + l)a 2 e- 1+ ~< 



+ - 



l n ^J + 0(l 



2c — v v 

,7 = 1 



dv + o(l) , 7 > 1 



7 < 1 



(61) 



(c - 6)^ +1 1 - 7 

If we are interested in quantities depending on the random potentials only via £, T and r 2 
as x(oo) and ^^(O), we can distinguish the following cases. 



1. 7 > 1. In this case r and r 2 have a finite limit as E j. i? c . As a consequence, i(c») 
is positive and Sa^O) is finite at E — E c and their value is given respectively by (pTf ) and 
(pl|): there is conduction with normal diffusion. So when E increases and goes through 
E c , both x(oo) and ^^(O) change discontinuously from to a positive value and then vary 
continuously. This is analogous with a first order phase transition. 

2. < 7 < 1. From the point of view of x(oo) the transition is still of first order, but 
the divergence of D, 



D 



E 3 £ 2 t 2 
2{i - A) 3 



(62) 



with the diverging r 2 when E \ E c resembles the divergence of the susceptibility at T c in 
second order magnetic phase transitions. Thus, in this case we have conduction accompanied 
with an anomalous diffusion. 



3. — 1 < 7 < 0. Both t and r 2 diverge when E j E c , so x(oo) tends to zero and 



D 



l 2 T 2 



2Y- 



oc e 



-1-27 



(63) 



with an additional factor | me| 3 if 7 = 0. So D diverges if 7 > —1/2, tends to zero if 
— 1<7<— 1/2 and to a finite nonzero limit if 7 = —1/2. 
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Let us emphasise that the probability density p 7 is purely continuous and, thus, the 
probability that fj = E c for a given j is zero. The probability that fj = E c for any j is still 
zero. So with probability 1 the particle will never be stopped, and x(oo) = means only 
that, with probability 1, x(t) = o(t), i.e. x(t) increases slower than t . It is in this way that 
we can understand the different possibilities of diffusion in the third case. 



V. SCALING AT CRITIC ALITY 

When E > E c , S&\&(u) is a meromorphic function of a; with no pole in a neighborhood of 
the origin. Analyticity at u — will be lost as E attains its critical value E c . In this section 
we derive a scaling law describing the behaviour of S±\±(u;E) when lj — ► and E { E c 
simultaneously in such way that z = e/oj = (E — E c )/oj is kept fixed. More specifically, we 
expect that in this limit 

S^u-E)^ 7 ^- (64) 
Below we prove the above form and find the scaling function r(z). The model that we use 



for explicit computations is the same as in section 4, given by ([55]) and fl5"7|), although the 
conclusions certainly hold more generally. First we would like to give an argument why one 
can expect the asymptotic form (1641). If e > 0, the correlation function 



£(t, if) = x(t)x{f) - x(t) x(t') (65) 

decays exponentially with \t — 1'\, and the correlation time can be approximated with r m = 
2aE/e(E + E c ), the maximum passage time through a scatterer (see equation (^6|) ) . This 
allows one to write down at least two different, qualitatively reasonable, approximations to 
£(t,f)> namely 

Ut,t') = ^(r m -\t-t'\) and f 2 (t, if) = £ e-^ T ™ (66) 
with Q(y) being the Heaviside function. Both lead to a form like (164]). From we obtain 



Sg(u; E) m — sinu;r m « ^ sin - (67) 



while ^2 yields 



q( 2 ) (uj . E ) ~ 2ioT ™ ~ 2e ° a,z (m 



In deriving r(z) we will consider only negative values of 7, so that < a = —7 < 1, 
choose 6 = for the sake of simplicity and use the notation 

P(y)= Py (E c -y)=py- a (69) 

cf. equation (0). We note that only 

sup y a P(y) < 00 and \imy a P(y) —p , (70) 
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where p is defined by the normalization, will be used below, so instead of (|69| ) we can take 
any P(y) with the properties (|7Dp. Notice that dy^—^- = oo. Accordingly, we obtain 



e c 



dyP(y) exp 



iuoa 



1 1 

+ 



e + y e + 2E c -y / 



(71) 



and 



lim ; = p I ds s 



exp % 



s + z 



1 



9{z) 



(72) 



Furthermore, from the third line of (|60|) we have 

t = r e- a + 0(e l - a ) (73) 
as e — > 0. In order to compute the current spectrum other terms have to be computed, like 

I _ e -/3M = ! _ e -i%('-A) ~ (74) 

Next, we need the average of 



C(/i) = /x r dye- 1 **® - £7(1 - e"" r ) 

J —a 



For the potential ( |55|) 
Inserting ( ff6|) in (|75|) , 



« ds 



CM = / 



1 — 2 exp 



/ia 



E + f , 



E + f ""E + f E-f 



1 1 
+ exp — /ia I — - + 



E + f E-f / 



so that 



coo 



+ exp 



dyP(y)(E c -y) l-2exp 
/ /ia /ia 



/ia 

'e + 2E c - V/ 



e + y e + 2E c - Vj 



and 



e^ r C(/i) 



dyP(y) (£ c -y) 
/ia 



/ /ia 
exp I — : h 



/ia 



—2 exp 

\ £ - v, 

Accordingly, in the limit of vanishing u 



+ 1 



(75) 



(76) 



(77) 



(78) 



(79) 
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0T 



roo 

^Ciiu) ^pE c dss- 
Jo 



exp l —i 



s + z 



and 



roo 

u a - x e^C(iu) ^pE c dss~ a 

Jo 

Moreover, one has (taking £ = a for simplicity) 



exp i 



s + z 



E c g*{z) 



E c g{z) 



B(-iu) 



C(-iuj) 



1 — exp (— (3(iu>)) 



uu a ~ l (l - exp (-/3(iu)) 
Similar expressions have to be obtained for 



- * 


E c g*(z)' 




g*(z) 



-E r . 



\C{iu)\ 2 = / dyP(y)(E c -yY 
Jo 



+ exp —i i- 



2 exp 



uoa 



exp —i 



e + y e + 2E c -y 
so that 

roo 

LO a ' l \C(iu)\ 2 ^pE 2 c / dss- a 

Jo 

By substituting into eq. ( |Hf ) we finally obtain 

e + E c 



V e + 2E C - y t 

2 



s + z 



- 1 



Ej.u(z) 



u 2 [£ - 2a + (e + E c )r] 
e + E c 



uj 2 uj a ~ l [t - 2a + (e + E c )t e 
Accordingly, in the limit of small e one obtains 



iuj)\ 2 + 2Re [e iujS B{^)e^C(hS) ] } 
E 2 C [u{z) - 2Reg(z)} 



E 2 



Lu a+L E- a S ±l± (iu) -> - 2Re^)] 

TO 



7? 2 „oo 

/ dtt~ 
T Jo 

4^2 



TO 



roo 

p / dtr a 

Jo 



1 - e"^ 
1 — cos 



2Re (e^ - 1 



t + z 



Finally, 



lim uSz\&(u; E c + ujz) 



AE 2 c z a J °°dtt- a 



1 — cos 



\t+z) 



a J oo dtt- a (l + t)- 1 



= r(z) 
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It is worth pointing out the relevance of this result: a uj^ 1 component emerges naturally at 
criticality in the spectral properties of a purely mechanical disordered model, that can be 
viewed as the classical counterpart of the Anderson's model for localization. 
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